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We discuss recent work on methods for incorporating nonlcading QCD corrections in 
parton shower algorithms. 



Parton shower Monte Carlo event generators are the main practical tool to de- 
scribe multi-particle final states in high energy collisions. These event generators 
couple a leading-order (LO) hard scattering to a showering, treated basically in the 
leading logarithm approximation. Many nonleading effects are also included in these 
event generators: for instance, through the exact multi-parton kinematics; through 
the angular ordering of gluon emission; through optimal choices of the renormal- 
ization scale 'in the running coupling. In some cases there are also treatments (see, 
e.g., refs.EJocra) to approximately include next-to-leading-order (NLO) corrections 
to the hard scattering. Nevertheless, there is as yet no method for going beyond 
the leading approximation systematically. This implies that event generators can- 
not incorporate fully the known NLO (and NNLO) calculations of hard-scattering 
cross sections. 

Recent papersBi'B have started to investigate systematic subtractive methods as 
a way to tackle this problem. This talk describes some of the work in this direction. 

To discuss this, let us schematically represent the cross section in an event 
generator as 

v[W]= W{X)P$®H. (1) 

final states X 

Here W is a weight function that specifies the definition of the particular cross 
section under consideration. The symbol PS denotes the parton shower and the 
symbol ® denotes its action on the initial and final partons in the hard scattering, 
whose cross section is denoted by H. 

In a standard Monte Carlo, the hard scattering is taken to the leading order, 
H( LO \ and PS denotes showering from the partons in H^°\ In an NLO Monte 
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Carlo, the cross section involves a structure of the form 



PS 



ff( LO > + a s ( ff( NL °) - PS/(1) ® ff (LO) - PSjp(I) ® ff( L °))l . (2) 



Here the first term in the square brackets is the LO hard-scattering function, and 
the second term is the subtracted NLO hard-scattering function. 7?( NLO ) is the 
result of computing the partonic cross section from the NLO graphs, while PS/(1) 
and PSf (1) are the order a s approximations to the initial-state and final-state 
showering. The subtraction terms avoid double counting of events already included 
by showering from H^ LO \ 

Ref.i gives a simple example of the structure (||): it treats the photon-gluon 
fusion process in deep inelastic ep scattering, using the Monte Carlo algorithm of 
ref.i. This case is special because i) only the subtraction term in PS/(1), coming 
from the initial state, is present, and ii) there are no leading-power contributions 
from soft gluons. Then the geometry of the leading-power regions that contribute 
to (ph is simple, just consisting of the ultraviolet region and the region collinear to 
the initial state, with no overlap. 

For general NLO processes, leading-power contributions to (||) involve a compli- 
cated geometry of possibly overlapping regions, that include soft contributions as 
well as collinear contributions. A typical case is shown in Fig. 1. Ref.Lj discusses 
how to extend the method to deal with such cases. A first step in this program is to 
develop techniques to decompose Feynman graphs into sums of terms over different 
regions, with the terms arranged so as to correspond to factors in a factorization 
formula suitable for the event-generator application. Similar problems are encoun- 
tered in multi-loop calculations based on graph- by-graph methods (see, e.g., ref.EI). 
This corresponds to a decomposition for i/( NLO ) of the kind 

^(nlo) _ ^ A H (R) + nonleading power, (3) 

regions R 

holding uniformly over the whole of the phase space. Each of the pieces in (|J) 
contains counterterms that prevent double counting and provide the suppression 
for going outside the region in which that particular piece was originally supposed 
to give a good approximation to the matrix element. This subtractive approach 
is to be contrasted with approaches based on splitting the phase space in different 
domains and using different approximations to the matrix element in these different 
domains (see, e.g., xefaJilH). 

Semi-analyticalEj'til or fully numerically subtraction methods have been devised 
to calculate NLO quantities that are infrared safe. These methods are not directly 
applicable in event generators that simulate the fully exclusive structure of the 
hadronic final states, since here the quantities being computed are not infrared safe 
in perturbation theory. In particular, one cannot use a cancellation of soft gluon 
contributions between real and virtual graphs. 

Ref.E3 discusses a strategy to construct a decomposition of the kind (||). This is 
inspired by the R-operation techniques of renormalization. See ref.Ej for a related 
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Figure 1: Geometry of leading-power regions for a typical NLO cross section. The 
axes denote lightcone directions in momentum space. 



approach. Given t 
counting arguments 



list of the leading regions, determined by standard power 
i3, we proceed from "smaller" to "larger" regions. See Fig. 2. 
For each region R, we remove the contribution from smaller regions, and construct 
an approximation to the matrix element valid in R, up to power suppressed cor- 
rections. Then we subtract any divergences that appear in this expression coming 
from larger regions. To ensure that the splitting between the terms is defined 
gauge-invariantly, at each step we demand that the counterterms be constructed 
from matrix elements of Wilson line operators, 

V{n) = "Pexp I ig / dyn ■ A(yn) I , (4) 



with suitable directions n for the lines. Evolution equations in n enable one to 
connect the results corresponding to different directions. Taking n along lightlike 
directions gives back the eikonals of the leading infrared approximation, while taking 
n off the light cone provides effective cut-off parameters — but in such a way that 
even a formalism involving off-shell partons (see, e.g., ref.u) could be treated with 
gauge invariance. It has beerishown how to apply this procedure to virtual loopalj 
and real emission corrections!! 

Ref.l3 uses this method to construct a decomposition (|^) for one-gluon emission 
graphs in deep inelastic scattering. The term corresponding to the ultraviolet region 
gives the subtracted hard-scattering function to be used in Eq. (|2|). The collinear 
terms correspond to the evolution kernels to be used in the showering. The soft term 
would correspond to a new element in the Monte Carlo, but refE shows that this 
term can be eliminated by a suitable choice of the directions n for the Wilson lines. 
This is a result analogous to one in ref.El Whether or not this result generalizes to 
all orders remains to be investigated and is very important to the construction of 
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4-dim. 




Figure 2: The soft (S), collinear (Ci, C2), ultraviolet (UV) regions, corresponding 
to the model geometry of Fig. 1. 



NLO parton showers. 

The decomposition of ref.LI entails a specific definition for the collinear factors, 
which will not necessarily coincide with the definitions used in any current event 
generator. The issue of deriving a showering algorithm that corresponds to the 
subtracted collinear terms is the subject of current work. A correct answer to this 
question will encompass soft-gluon coherence and angular ordering. It will involve 
evolution equations with respect to the direction of the Wilson line, Eq. (Q), in 
terms of which the collinear subtractions are defined. 
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